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ABSTRACT 

In Part I of this paper, an overview of parabolic and PNS 
methodology developed to treat highly curved sub and supersonic wall 
jets is presented. The fundamental data base to which these models 
have been applied is discussed in detail. The analysis of strong 
curvature effects has been found to require a semi-elliptic exte nsio n 
of the parabolic modeling to account for turbulent contributions (v'v’) 
to the normal pressure variation, as well as an extension to the 
turbulence models utilized, to account for the highly enhanced mixing 
rates observed in situations with large convex curvature. A non- 
iterative, pressure-split procedure is shown to extend parabolic models 
to account for such normal pressure variations in an efficient manner, 
requiring minimal additional run time over a standard parabolic 
approach. Curvature corrections to a ke two-equation turbulence model 
are reviewed and their general applicability is assessed. For complex 
flows, the use of algebraic or full Reynolds stress turbulence models 
may be required, but the ke corrections utilized have been adequate for 
all fundamental cases thus far explored. For strong blowing situa- 
tions, a supersonic/underexpanded wall jet structure develops with a 
complex multiple shock cell internal wave structure. A new PNS 
approach is presented to solve this problem which extends parabolic 
methodology via the addition of a characteristic-based wave solver. 
Applications of this approach to analyze the interaction of wave and 
turbulence processes in wall jets are presented. The present uncer- 
tainty in dealing with compressibility effects in supersonic problems 
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is pointed out as a problem area for which no data exists. The 
unification of the parabolic, pressure-split and PNS wave solving 
capabilities into the wall jet computer code, VJET, is discussed. This 
code has served as a research tool for studying the effects of various 
parameters on wall jet structure, and includes advanced turbulence 
models with curvature and compressibility effects. In Part II of this 
paper, the steps taken towards incorporating WJET into a zonal 
component model for analyzing circulation control airfoils is 
presented. 


IHTBODUCTIOH 

A zonal model (Figure 1) for the engineering analysis of circula- 
tion control airfoil performance (TRACON) was developed by Dvorak and 
coworkers under David Taylor R$D Center (DTRDC) support 1 ,*. TRACON is 
comprised of separate components which analyze the external potential 
flow (Jameson's FL06 model), the airfoil boundary layer (Cohen and 
Reshotko /Green, laminar/ turbulent integral models), and the wall jet 
(Dvorak's finite difference model*). 


EXTERNAL (POTENTIAL) 

BOUNOARV LAYERS FLOW STREAMLINES 



FIGOBE 1. Zonal Approach for CC Airfoil Analysis. 


While TRACON was demonstrated to perform quite well for a variety 
of cases, limitations in its ability to deal with very large curvature 
and/or strong blowing were encountered. These limitations were asso- 
ciated with the modeling assumptions in the TRACON wall jet component 
which include: 


(1) a parabolic approach does not solve the normal momentum 
equation across the viscous wall jet. hence - AC is not 
predicted and has to be estimated from inviseid 
considerations; 

(2) a parabolic approach which does not treat the supersonic 
wave/shock structure occurring in underexpanded wall 
jets at high rates of blowing; 

(3) a weakly interactive, displacement-thickness based 
viscous/inviscid coupling approach which becomes in- 
adequate for strongly interactive situations associated 
with large curvature and/or strong blowing; and, 

(4) an algebraic eddy viscosity turbulence model which does 
not handle lag effects associated with significant 
pressure gradients. 


To remedy these limitations, a new wall jet model, WJET, was 
developed by Dash and coworkers (under DTRDC support) , which employed 
advanced numerical procedures and utilized a two-equation turbulence 
model. The first version of this model 4 ,* solved the higher order 
parabolic curved wall jet equations utilizing a conventional implicit 
algorithm. Ihis model provided for a formal solution of the viscous 
normal momentum equation to yield ACp across the jet, and, employed a 
hybrid (inner VanDriest/outer ke) two-equation turbulence model with 
curvature correction terms. An improved version of WJET*, 7 utilized a 
subsonic pressure— split approach which solved the semi-elliptic coupled 
continuity and normal momentum equations across the jet in a non- 
iterative manner, and, provided for direct coupling of the wall jet 
with a potential external flow solver, eliminating the requirement for 
weakly interactive displacement thickness based coupling. The final 
version of WJET developed under DTRDC support*,*, extended the 
parabolic/pressure-split methodology to provide PNS spatial marching 
capabilities in supersonic flow regions. A new implicit/explicit 
approach was utilized which employs an upwind finite difference 
representation of viscous-characteristic methodology to solve the wall 
jet wave field. 

In concurrent work initiated in 1984 under NASA/Ames support, the 
methodology was formulated 10 and is now being made operational, to 
incorporate WJET into the TRACON code, replacing the existing wall jet 
component. This effort is being performed jointly by Dvorak and 
coworkers (at AMI) and Dash and coworkers (at SAIC) . Progress towards 
this end is discussed in Part II of this paper 11 . 

In Part I of this paper, a brief overview of the features and 
capabilities of WJET, and, it 3 application to various simple cases will 
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be discussed. WJET has been used on a stand alone basis to analyze a 
variety of basic wall jet data which has led to a greater understanding 
of the dominant Influence of turbulence modeling on our ability to 
simulate wall jet flowfields. On highly curved convex surfaces, con- 
ventional two-equation models grossly underestimate the Jet growth and 
the ng that occurs, and curvature correction terms are required to 
properly simulate the flowfield. At supersonic velocities, free jet 
data indicate that compressibility effects can markedly reduce Jet 
growth and mixing, as will be discussed. A good data base to isolate 
the influence of compressibility effects on wall jets is not available, 
which leads to some uncertainty in our ability to treat circulation 
control airfoil problems, but a much larger uncertainty in other higher 
speed problems (i.e., tangential injection problems in supersonic 
combustors, slot cooling problems for hypersonic vehicles, etc.). 


OVERVIEW OF WJET MODEL 
Mean Flow and Turbulence Model Equations 

WJET solves the higher order curved boundary layer equations listed 
in Table 1. The equations are east in surface-oriented s,n coordinates 
(Figure 2) and include a tracer species equation for 6 (d = 1 in 
unmixed jet, = 0 in airstream) to delineate the jet/air mixing region. 
A classical Boussinesq approximation is utilized to relate turbulent 
stress terms to mean flow gradients, with the parabolized stress terms 
retained listed in Table 2. Turbulence closure is achieved using the 
two— equation ke model with standard coefficients (C = 0.09, C x — 1.43, 
C a = 1.92, a = 1.0, a = 1.3). The turbulent transport equations for 
k and e are Tisted in &ble 3. 

To extend the high Reynolds number ke turbulence model to the wall, 
a variety of techniques are available ranging from simple wall function 
approaches to the use of low Re extensions of the ke model. A review 
of these techniques from both a pragmatic and computational viewpoint 
(see ref. 5) has led to our use of a classical VanDriest mixing length 
formulation with damped law of the wall, as routinely employed in two 
layer algebraic turbulent model formulations (i.e., this comprises the 
inner layer component in the popular Cebeci— Smith and Baldwin— Lomax two 
layer formulations). Coupling between the inner (near wall) mixing 
length formulation and the outer ke formulation is set to occur at 
y = so. The values of k and e at the matching point are determined 
via the requirement that the mixing length and ke turbulent viscosities 
match and that the turbulence is in a state of equilibrium. This 
yields lower boundary conditions for k and e at the matching point. An 
analogous ML/ke coupling procedure has been developed by Arora et.al 1 * 
for application to a variety of turbulent boundary layer problems. The 
inner /outer coupling relations are listed in Table 3 along with a 
schematic of the coupling procedure. Results obtained are relatively 
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TtRT.it 2 - Parabolised Stress Terms 


Using a Boussinesq type approximation 
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Insensitive to the value of y + utilized for the ^oupling as long as it 
nominally remains in the log region (viz., 20 ~ y ~ 100). 



K- 1 /* K + 1 /* 

SURFACE-ORIENTED COORDINATES MAPPED GRID NETWORK 


FIGURE 2. Surface-Oriented Grid Nomenclature and Mapped/Stretched 

Grid Utilized. 


While both the mean flow and turbulence model equations contain a 
number of curvature terms arising from the transformation to curvi- 
linear, surface-oriented coordinates, numerous investigators have 
demonstrated that additional, curvature corrections terms are required 
to account for the strong effects of curvature on wall jet turbulence 
structure. The analogy drawn by Bradshaw 1 * between curvature and 
buoyancy has been utilized by most investigators as the basis for 
heuristic corrections to algebraic or two-equation turbulence models. 
Defining the curvature parameter, s = -KU/OU/3n), a curvature correc- 
tion to the ML formulation is given by: 
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where L is the planar length scale value, while a is a constant, (5 < 
a < 10)°. Ibis treatment has been implemented for curved wall jets by 
Foloyan and Whitelaw 14 who utilized a complete (inner/outer) mixing 
length formulation. For the present near wall use of the mixing length 
formulation, this correction will only be required in situations with 
very large curvature. For problems with small to moderate curvature, 
the near wall region correction to L q is negligible. 

Launder and coworkers 15 have developed a curvature correction for 
the ke model which utilizes a single empirical coefficient, C^. The 
curvature correction is proportional to a Richardson number, Ri, based 
on the turbulence time scale. In their formulation, the local Richard- 
son number is given by: 


Ri = -KO ( — )* {a 
6 on 


( 2 ) 


and the C a coefficient of the e equation is modified as follows: 


C, = 1.92 (1 - C Ri) 

2 c 


(3) 


Values of C of about 0.2 have yielded optimal predictions for a 
variety of cuWed boundary layer flows as described in reference 15. 

An analogous type of curvature correction for the ke model has 
been developed by Hah and Lakshminarayana at Penn State 14 . They have 
modified the C x coefficient of the e equation as follows: 

C, = 1.43 (1 + C Ri) (4) 

1 c 


We have implemented the Penn State correction using the Launder defi- 
nition of Ri in our model. 

At supersonic wall jet velocities, a compressibility correction 
analogous to that of Dash, et.al. 17 for free shear layers, may be 
required. Here, a factor whose magnitude varies from 1.00 at M * 1 to 
0.25 for large Mach numbers is used to multiply the C coefficient of 
the ke model. This correction factor is an empirical expression derived 
from experimental free shear layer observations and is calculated as 
follows : 
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( 5 ) 


f = 0.25 + 0.75/ {1.0+exp[24.73 (M -0.2) ] J 

cc t 

where M is the Mach number characterizing the fluctuating velocity 
field, fe.g., k 1 /* divided by the local speed of sound). 

For free shear layers, M is evaluated at the position of maximum 
k at each station. Its adaptation to supersonic wall jets would entail 
utilizing the value of k at the position of maximum velocity and 
incorporating a transverse dependence to smoothly reduce the correction 
in the 'boundary layer' portion of the wall jet (below the peak 
velocity point). Data exhibiting the dependence of mixing on the wall 
jet Mach number is not presently available to derive a correlation 
analogous to that of equation (5) for wall jet flows. 


HUMEBICAL PROCEDURES 

Splitting of Solution into Parabolic and Elliptic/Hyperbolic Components 
The approach taken in WJET involves combining: 

(1) a parabolic solution of the streamwi3e momentum, energy, 

species parameter, and turbulence model equations with 
the streamwise pressure gradient term, 3P/3s (s,n) 

specified - this solution yields the variation of U, H, 
6 , k and e; 

(2) an elliptic/pressure-split solution of the coupled 
continuity and normal momentum equations in subsonic 
regions which yields the variation of pressure and 
normal velocity across the wall jet; and, 

(3) a hyperbolic/upwind characteristic-based solution of the 
coupled continuity and normal momentum equations in 
supersonic regions which yields the local pressure and 
flow angle. 

These three solution procedures are unified in the WJET code to provide 
generalized spatial marching capabilities for a broad category of wall 
jet problems. 

Parabolic Analysis 

The WJET parabolic algorithm integrates the U momentum, H, 6 , k 
and e equations (Table 1) in mapped rectangular coordinates (Figure 2). 
The mapped, vectorized equations take the form: 
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( 6 ) 


where: 


f - [ 0 , H, 6 , k and e ] T 

a and b are mapping parameters, and g f is the source term. The equa- 
tions are spatially integrated using an upwind/implicit algorithm. A 
fixed number of grid points are distributed between the wall (n ■ 0) 
and outer viscous boundary ( t| * 1) whose growth is obtained via adap- 
tive methodology keyed to the edge gradients. The distribution of grid 
points, t| ( I ) , remains invariant throughout the calculation and the 
stretching utilized can be arbitrarily stipulated, or specified using 
built in grid distribution parameters. The equations are solved in an 
uncoupled manner (the source terms are solved explicitly) and the 
difference equations then take standard tridiagonal form. Complete 
details of the parabolic algorithm are available in references 4 and 5. 

Pressure-Split Subsonic Cross-Flow Analysis 

To analyze subsonic wall jets with large curvature, a pressure- 
splitting approach analogous to that of Bradshaw and coworkers 1 * , ±9 is 
utilized. In the pressure-splitting procedure, the global pressure 
field, P*(s,n), utilized to evaluate dP/0s in the streamwise momentum 
equation, must be stipulated. This is initially estimated to be the 
inviscid pressure field prevailing in the region occupied by the wall 
jet. In the pressure-splitting approximation, the parabolic equations 
are integrated with 3P/ds obtained from P*(s,n). However, the pressure 
field is revised in the course of the spatial integration by solving 
the coupled continuity and normal momentum equations across the jet 
with the inviscid pressure prevailing at the edge of the jet serving as 
an outer boundary condition. Global convergence in regions of strong 
curvature is obtained by repeating the calculation with the revised 
pressure field until the imposed and upgraded pressures are effectively 
the same. 

By manipulations described in references 6 and 7, the continuity 
equation (in mapped coordinates) can be written in terms of pressure, 
P, and normal velocity, v, yielding: 



( 7 ) 
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where the source term. 



is listed below: 
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Note that g contains partial derivatives of U and H (which are 
evaluated a jlriori as part of the parabolic solution) and also contains 
the prescribed streamwise pressure gradient term. The normal momentum 
equation given below (g represents turbulent stress terms - see 
Table 2): Y 




+ bh + 


K pU = g 


( 8 ) 


is solved with the continuity equation in a coupled manner to yield the 
variation of P and v across the wall jet (see references 6 and 7 for 
details) . 

Figure 3 illustrates results obtained using this pressure-split 
approach to analyze the simple case of a curved wall jet issuing into 
still air. Shown are the wall jet geometry, maximum velocity decay 
(contrasted to the variation for a flat wall), induced entrainment 
(also contrasted to the flat wall variation), surface pressure vari- 
ation (utilizing viscous and inviscid forms of the normal momentum 
equation) and a decomposition of terms in the normal momentum equations 
showing their individual contributions to AC across the wall jet. 
Note that the very significant contribution of turbulent stress terms 
(specifically v'v') to the AC across the jet. Complete details of 
this case are provided in references 6 and 7 . 

Upwind Characteristic-Based Supersonic Wave Solver Analysis 

To analyze supersonic regions of underexpanded wall jets, a 
characteristic-based procedure is employed to locally evaluate the 
wave field (e.g., to obtain pressures and flow angles at each grid 
point). The approach taken involves a 'modern' formulation of viscous- 
characteristic methodology originally developed about 20 years 
agoio^i. The viscous-characteristics approach involves manipulating 
the continuity and normal momentum equations to obtain characteristic 
relations of the form: 


3 1 * 
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FIGURE 3. Pressure-Split Analysis of Curved Wall Jet Issuing into 

Still Air. 


2iae£22lt ± i9 . -Jtlrt.tnf , p± dl ± (s) 

where the source term, F, contains the viscous stress terms, and 
diffusive transport terms, appearing on the r.h.s. of the streamwise 
momentum (F normal momentum (F ) , and energy (F.,) equations . The 
source term, F, is given by: u 


F = A(F u +F y ) + BF r (10) 


where : 

A = -[(l+(y-l)M*)sinjicos9 ± cos|isin0]/ (yPM* ) 

and: 

B = -(y-Dsinp/ (yPQ) 

The nomenclature for the characteristic equations is exhibited in 
Figure 4. 



FIGURE 4. Nomenclature for Characteristic Equations. 
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In earlier applications of this approach, inverse characteristic 
methodology was utilized which was extremely cumbersome. The new 
approach developed*,* involves representing total derivatives along 
characteristics as a combination of streamwise and normal partial 
derivatives in a mapped (5,i\) computational space. The normal deriva- 
tives are evaluated at the known station using an upwind formulation 
keyed to the X - characteristic direction. With the manipulations 
described in references 8 and 9, the pressure at grid point I at { + A( 
can be evaluated as a function of the pressures and flow angles at grid 
points 1-1, I and 1+1 at station 5* This new formulation is summarized 
in Table 4. In present applications of WJET, the coupled parabolic/ 
hyperbolic solution is performed in the following three step sequence: 

(1) Prediction of wave field (P,0) at 5 + A? solving 

pressure equation and comparable flow angle equation 
(Table 4) using coefficients evaluated at § at charac- 
teristic intersection points, and, viscous terms, F , F y 
and Fg evaluated at $ at grid points I. 

(2) Solution of parabolic system of equations (eq. 6) 

yielding f at £ + A2j. Pressure gradients 'prescribed' 
in accordance with wave field solution of step (1) . 

(3) Correction of wave field at § + using coefficients 

averaged along characteristics and, values of F , F and 
F„ which are evaluated at 5 + if the parabolic 

algorithm of Step (2) is fully implicit; or, are 

averaged across the integration step if a Crank-Nicolson 
parabolic procedure is utilized. 

The formulation was first checked out in the inviscid limit for 
weak shock-capturing capabilities and produced results comparable to 
those of the SCIPPY code (explicit MacCormack algorithm) as exhibited 
in Figure 5 (see ref. 8 for details) . For strong shock waves in 
inviscid regions, artificial viscosity must be introduced to stabilize 
the shock calculation and to generate entropy. The approach taken 
parallels that utilized in Beam and Warming based PNS algorithms and is 
described in reference 22. 

The wave solver formulation was then checked out for supersonic 
viscous /inviscid jet interaction problems by comparing WJET predictions 
with those of the well tested SCIPVIS PNS jet mixing model (see refs. 
23 - 26). These comparisons are described in references 8, 22 and 27. 
A typical comparison is exhibited in Figure 6 showing the interaction 
of an expansion fan with a free turbulent shear layer. 












Unified Parabolic/Preaaure— Split/Wave Solver Model 

The WJET code unifies the three solution procedures described 
above, providing for PNS-based spatial marching capabilities comparable 
to those developed for free jet problems (see Table 5 and refs, 23 — 
26). Referring to Figure 7, the wall jet problem is more complicated 
than the corresponding free jet problem since subsonic/supersonic 
coupling is required at both inner (I*) and outer (1$) sonic lines. 



FIGURE 7. Flow Segmentation for Underexpanded Wall Jet Problem. 


For free jets, a data base exists for underexpanded problems (see, 
e.g., ref. 28) which has been utilized for detailed verifica- 
tion**, **,*‘ . Figure 8 exhibits predicted wave/mixing layer structure 
of mildly underexpanded Mach 2 free jet issuing into still air, and, 
c ompa risons of predicted and measured turbulent Intensities, u'u' 
(u'u' = 2gk where k is the predicted turbu lent kinetic energy and g is 
an isotropy parameter used to related k to u'u'; g = 2/3 represents the 
isotropic situation and appears to best correlate with the measure- 
ments). Figure 9 exhibits comparisons of the predicted pressure 
variations (axis and off axis) with the data. The comparisons exhibited 
here (and the additional comparisons described in refs. 23, 25, 26 and 
28) are quite good and were obtained using the kW turbulence model. 
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TABLE S - Aspects of PMS Approach for Free Jets 



• Shock-Capturing PNS Solution in Supersonic Mixing 
Regions 

• Pressure-Split Approximation in Subsonic Mixing Regions 

0 Subsonic/Supersonic Coupling at Viscous Sonic Lines 

• Dse of Compressibility Corrected Two-Equation Turbulence 
Models 

• Solution in Mapped Coordinates Encompassing Viscous/ 
Inviscid Jet 

0 Direct-Coupling with External Flow Solution at Outer 

Viscous Boundary 
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MACH 2 JET; SCIPVIS, KW PREDICTION 




FIGURE 8. Predicted Wave/Mixing Layer Structure of Underexpanded 

Mach 2 Free Jet and Comparisons of Turbulent Intensity 
with Data. 
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PRESSURE VARIATION 



FIGURE 9. Comparison of Predicted and Measured Streamwise Pressure 

Variations; kW Model Solution. 
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Results with a standard ke model exhibit too fast a rate of mixing as 
would be expected from earlier comparisons with simple (balanced 
pressure) shear layer and jet data at supersonic velocities. The issue 
of compressibility effects in supersonic free turbulent shear layers 
has been addressed from a pragmatic viewpoint (see refs. 17, 29 - 31) 
and two-equation turbulence models are now available’* which can 
adequately analyze the rather broad base of 'fundamental' high speed 
free shear layer/ Jet data and also, some non-fundamental situations with 
significant wave structure, as exhibited above. 

Preliminary calculations made with WJET for an underexpanded 
curved wall jet are exhibited in Figures 10 through 12. Figure 10 
exhibits the overall geometry and predicted streamwise variations of 
principal jet surfaces (viz., inner /outer shear layer boundaries where 
4 - 0,95/0.05, jet half radius, outer sonic line, and outer/adaptive 
computational boundary) in physical (x,y) and surface-oriented coordi- 
nates. Figure 11 exhibits the variation of wall pressure and skin 
friction coefficient - note the rapid response of skin friction to the 
wave field (the details of the interactive procedure for analyzing the 
near wall subsonic portion of the wall jet are described in ref 3 . 8 and 
9) . Figure 12 exhibits predicted profiles of pressure and Mach number 
across the jet. Also shown is the normal grid distribution which is 
highly stretched in the near wall + region (viz., from wall to position 
of velocity maximum - typically, y of second grid point is ~ 1 and the 
same number of grid points span the near wall region and the outer 
region), and, equally spaced in the outer region (from maximum velocity 
position to outer boundary) . 

Unfortunately, adequate data to validate underexpanded wall jet 
solutions and thus resolve turbulence issues regarding compressibility 
effects, etc., is not presently available, and, hence, no such compari- 
sons with data are exhibited in this article. The high speed wall jet 
data base available has recently been reviewed (under programs geared 
towards high speed film cooling and tangential injection in supersonic 
combustors) . No data has been identified as suitable for turbulence 
model validation due to lack of key measurements (e.g., detailed initial 
profiles, turbulence levels, etc.). These issues are described in 
references 33 and 34. 

Coupling Procedures for Hall Jet and External Potential Flow 

To incorporate the wall jet code in a zonal approach for analyzing 
circulation control airfoils, a variety of coupling techniques were 
reviewed as discussed in references 4 and 5. Figure 13 schematizes 
coupling procedures available for subsonic wall jet 3 . The displacement- 
thickness coupling approach overlaps the inviscid and viscous solutions 
and utilizes standard boundary layer concepts which break down for 
thick, highly curved viscous layers where the normal pressure variation 
is significant. The direct pressure-split coupling approach introduced 
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FIGURE 10 


Predicted Features of Curved Underexpanded Wall Jet 
Issuing into Still Air. 
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FIGURE 12. Pressure and Mach Number Profiles at Position Fifty Slot 
Heights Downstream of Slot. 


by Bradshaw and coworkers 1 •, 19 couples the viscous and inviscid solu- 
tions at the jet viscous boundary, and the pressure field within the 
jet is determined by the pressure-splitting methodology described 
above which directly accounts for the contribution of stress/diffusive 
terms. Details of this coupling methodology for wall jets are 
described in references 6 and 7. Velocity-split coupling also directly 
involves a complete overlap of viscous and inviscid solutions, but 
here, the coupling is intimate and can account for separated flow 
regions. Applications of this approach to nozzle afterbody problems 
have been quite successful (see refs. 37 and 38) and results comparable 
to full NS results have been achieved in a fraction of the run time. 



DISPLACEMENT- THICKNESS COUPLING 


COUPLING 


WALL JET 


POTENTIAL FLOW 


POTENTIAL flow 





WALL JET 


V e along n # — £ n along n s b 

J^.Ue.Tfl - P b .U b ,T b - 
0 s 


DIRECT PRESSURE- SPLIT COUPLING 

I 


POTENTIAL FLOW 


WALL JET 


V e = *n 


COUPLING 


• As Above but Match along ne 

• P(s, n) in Jet from Normal - 
Momentum via Several Iterative 
Sweeps for Given Pe 


VELOCITY-SPLIT COUPLING 


POTENTIAL FLOW 


♦ U = u(s, n) 

WALL JET 


COUPLING 
U/U E = u ( I + <£ s ) 

V/Ue = <£n 

• Solutions Overlap 

• Solve Equation with u(s,n) Prescribed 

• Solve Wall Jet Equations with <f> 
Prescribed 


FIGDHE 13. Coupling Procedures for Subsonic Wall Jet. 


For supersonic wall Jets, the viscous/inviscid coupling require- 
ments become more complex and no coupling methodology has as yet been 
made operational. For supersonic free jets, the overlaid viscous/ 
lnvlscid coupling approach of Dash, et.al. 10 , which employs displace- 
ment-thickness coupling concepts 40 (Figure 14) has been successful In 
studies geared toward nozzle afterbody drag predictions 41 . The RAXJET 
zonal component model (which is based on this methodology) employs 
components (Figure 15) analogous to those utilized In the TRACON CC 
airfoil code. The extension of the overlaid coupling approach to wall 
jets was investigated by Dash 41 but found not to be a viable method 
(see references 4 and 5). 

Recently, the direct-coupling approach for supersonic free Jets 
has been made operational 41 utilizing a free jet version of WJET (the 
SPLIT? model 17 ) coupled to the VSAERO panel method potential solver 44 . 
The work (supported by AFWAL) is geared towards developing interactive 
methodology for VSTOL Jets (Figure 16). Typical predictions are 
exhibited in Figures 17 and 18. Figure 17 depicts the direct-coupling 
boundary comprised of the nacelle surface and paneled jet boundary, 
y B , which lies close to the outer jet computational boundary. Also 
exhibited is the source distribution, 6 , applied along y B , which 
combines the effects of jet entrainment n ( suction) and jet Dlockage 
(underexpansion - shock effects) . Figure 18 ^depicts the predicted 
pressure variation along the coupling boundary y fi . In the first pass 
iteration, the jet is represented as a solid sting and the pressure is 
given by C( l ). In subsequent iterations, the coupled effects of jet 
entrainment? and blockage are evident in the pressure variation. See 
references 29, 43 and 45 for further details of this methodology. 
Progress towards using advanced coupling concepts for incorporating 
WJET into the TRACON CC code 11 , is described in Part II of this 
paper 11 . 


AHALTSIS OF FUNDAMENTAL DATA 

The data analyzed by WJET are limited to situations for which WJET 
can operate on a stand alone basis. Hence, the cases involve rather 
fundamental situations and primarily reflect upon the ability of the 
turbulence model incorporated in WJET. The analysis of wall jets in 
realistic CC airfoil flowfields requires coupling of WJET with 
TRACON 11 . Most of the cases analyzed have already been described in 
refs. 4 and 5 and only a very brief overview will be provided here. 

Planar Wall Jet Issuing Into Still Air 

For this simplest of all wall jet cases, the overall jet growth 
parameters (viz., half radius and locus of maximum velocity) predicted 
using the hybrid ke/VanDriest turbulence model are in reasonable 
agreement with the data (see Figure 19), and, in better agreement than 
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x Subsonic/transonic tilernol flow map 
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FiGUHE 14. Overlaid Viscous/Inviscid Coupling Approach and 
Displacement- Thickness Representation of Plume Boundary. 
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RAX JET COMPUTER MODEL (WILMOTH, NASA TM 83235, 1582) 

South/ Jameson Potential Flow Solver 
Integral BL Mooel 

SCIPPY Shock -Capturing Plume Mooel (Dash/Thorpe) 

. BOAT Overlaid Shear Layer Solution (Dash/Pergament) 

FIGURE 15. Components Utilized in RAXJET Code for Zonal Solutions 
of Nozzle Afterbody Drag. 


51 






52 



FIGUBE 18. Coupling Boundary Pressure Variation for Several Shock 
Cells . 


that of other investigators (e.g., see ref. 46) who used the ke model 
with a wall function near wall approximation. Also exhibited in Figure 
19 is a comparison of predicted and measured maximum velocity decay 
with the data correlation of Rajaratnam 4 7 . Figure 20 compares 
predicted and measured 4 * streamwise and normal velocity profiles at 
selected stations downstream of where similarity is achieved. 

Planar Vail Jet with Moving Outer Stream 

A number of calculations were performed 4 , 7 corresponding to 
experiments performed by Kacker and Whitelaw 4 *. Figure 21 exhibits 
typical comparisons achieved for maximum velocity decay, half radius 
and maximum velocity locus variation, and wall skin friction variation. 
The comparisons are quite good. 
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Cured Vail Jet lagging Into Still Air 

Calculations were performed',* corresponding to the experiments of 
Wilson and Goldstein 4 *. Figure 22 depicts comparisons of half radius 
variation and maximum velocity decay utilizing the standard (curvi- 
linear, s,n) ke model and versions with curvature correction 
terms 1 *, 1 *. The standard ks-based prediction is seen to grossly 
underestimate the rate of mixing while curvature correction predic- 
tions (using the recommended curvature coefficient-based values of C = 
0.2 [Launder 1 *] and C q = 0.16 [Penn State 1 *]) agree quite well with ?he 
data. Profiles of turbulent shear stress performed with the two 
correction terms (Figure 23) agree reasonably well with each other and 
with the data at 0 = 90°, but diverge at 0 = 180°, a 3 exhibited. 
Figure 24 exhibits the streamwise variation of peak turbulent shear 
stress and clearly exhibits the divergence in the predictions at 0 = 
90°. The data supports the Penn State correction, except for the 
abrupt jump at 0 - 180°. 

Carved Vail Jet Vith Moving Outer Stream 

The last calculation simulates one of the experiments of Kind* 0 , 
as schematized in Figure 25 - conditions correspond to the Flow II Case 
listed. The calculation was run using the pressure-split approach with 
conditions (pressure, streamwise velocity) prescribed at the jet outer 
edge. The predicted AC across the jet (Figure 26) is in very good 
agreement with the data except for 0 > 60°. The global pressure 
iteration approach was employed to eliminate the pressure-split 
approximation and after several iterative sweeps, significant improve- 
ment in the comparisons was obtained. The predicted variation in 
maximum velocity decay is exhibited in Figure 27 and the results with 
the curvature correction are significantly better that those with the 
standard ke model. The predicted variation in jet half radius is 
exhibited in Figure 28 and again, the improvements utilizing the 
curvature correction are quite significant. 


COHCLUDIHG REMARKS 

The ability to analyze fundamental wall jet data is clearly keyed 
to the capabilities of the turbulence model utilized. Our starting 
point had involved the use of a high Reynolds number two-equation ke 
model with heuristic corrections for streamwise curvature. This model 
was coupled to an inner (near wall) damped VanDriest mixing length 
model. Algebraic (eddy viscosity) models were not utilized since they 
cannot readily deal with complex jet/boundary layer length scales, with 
initial (slot /boundary layer) turbulence levels, and with lag effects. 
The curvature modifications to the e equation 1 *, 1 *, previously demon- 
strated to yield improvements for curved boundary layers, also appear 
to work quite well for curved wall jets. The use of these curvature 


57 


20 

16 

12 

y '/,/b 

8 

4 

Q 


L0 

.8 

JS 





FIGORK 23. Profiles of Turbulent Shear Stress for Curved Wall Jet 
into Still Air. 
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corrected models is recommended as a logical starting point for 
inclusion in CC NS solvers, since their ability to analyze fundamental 
wall jet data is reasonably well established. For supersonic wall 
jets, some type of compressibility correction may be required to deal 
with the near slot shear layer effects, but data is not presently 
available to support the heuristic modeling of such a correction. 
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Prediction of AC Across Wall Jet for Kind CC Airfoil 
Problem. ^ 
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